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SUMMARY 


This work is concerned with the self-consistent calculation of the electronic energy 
levels of noble gas pseudomolecules formed when a metal surface is bombarded by 
noble gas ions and with the construction of energy level correlation diagrams as a func- 
tion of interatomic spacing. The self-consistent-field xa scattered wave (SCF-Xa-Sl^ 
method of Slater and Johnson is utilized. Preliminary results on the Ne-Mg system are 
given. 

An interactive xa programming system, which has been implemented on the LeRC 
IBM 370 computer, is described in detail. This automated system makes use of special 
PROCDEFS (procedure definitions) to minimize the data to be entered manually at a 
remote terminal. Listings of the special PROCDEFS and of typical ii^ut data are 
given. 


INTRODUCTION 

Ferrante and Pepper (ref. 1) recently investigated Auger electron production by 
Ne'*’ bombardment of Mg, Al, and Si surfaces. Ion energies in the range of 
1-3 KeV were used. Auger electrons characteristic of both the target atoms and the 
incident Ne were observed. They attributed these observations to asymmetric 
Ne-surface atom collisions, and presented evidence that the Auger electrons resulted 
from autoionization of neutral, excited Ne. Their results on the variation in strength 
of Ne Auger signals with target in the sequence Mg, Al, Si appear to be in agree- 
ment with the predictions of the electron-promotion excitation mechanism of Barat and 
Lichten (ref. 2). Ferrante and Pepper (ref. 3) have also investigated the production of 
Auger electrons by bombardment of Mg and Al surfaces with other noble gas ions. 
To help elucidate the question of electron-promotion in the formation of collision- 
induced pseudomolecules, it is of great interest to obtain accurate energy level correl- 
ation diagrams for the pseudomolecules formed from various noble gas-metal pairs. 

In the present paper the results of a preliminary calculation of a Ne-Mg correla- 
tion diagram using the xof method are presented. In addition, a description of the 
xo program, which can be used for a wide variety of chemical and solid state calcula- 
tion, is presented to facilitate its use by interested parties at LeRC. 

Professor J. G. Norman, Jr. , of the University of Washington, supplied basic xa 
programs and operating instructions. 
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SELF-CONSISTENT CALCULATION OF PSEUDOMOLECULE ENERGY 
LEVELS, CONSTRUCTION OF CORRELATION DUGRAMS 

This work is concerned with the self-consistent calculation of the energy levels of 
noble gas- metal pseudomolecules formed when a metal surface is bombarded by noble 
gas ions, and with the construction of energy level correlation diagrams. The theoret- 
ical method utilized in the molecular calculations is the self- consistent- field xa 
scattered-wave (SCF-xa-SW) method of Slater and Johnson (ref. 4). The interactive 
xQf programing system described in the following section was in^)lemented on the 
LeRC IBM 370-3033 computer for this purpose. To construct a correlation diagram 
for a particular system, electronic energy levels must be calculated for a number of 
different interatomic spacings ranging from the separated atom to the united atom 
limits. The atomic calculations were performed using the modified Herman Skillman 
SCF program described below (ref. 5). Preliminary results on the Ne- Mg system are 
summarized in figure 1. 

In the usual application of the xof-SW method, the molecule is surrounded by a 
sphere, and the interior is partitioned into a number of spherical atomic regions and 
an interatomic region. The molecular potential is constructed to be spherically sym- 
metric in the atomic and extra molecular regions and constant in the interatomic 
region. The Schrodinger equation is then solved self- consistently, subject to the 
boundary condition that the wave function and its first derivative are continuous across 
all boundaries, using the multiple scattering formalism combined with the xa density 
function formalism for treating the exchange and correlation. In general, the iterative 
self-consistent field computational scheme is reasonably rapidly convergent. In the 
earliest applications of the xa-MS method, the atomic spheres were taken as non- 
overlapping. As larger and more "open” systems were investigated, the constant 
potential interatomic regions become an appreciable portion of the molecular volume, 
and convergence difficulties were found. This was remedied empirically by allowing 
the atomic spheres to overlap each other. Consequently, the theory was reformulated 
for the case of overlapping spheres. The major constraint on the use of overlapping 
atomic spheres is that the spheres may not overlap any adjacent atomic center. Thus, 
the largest allowable sphere for a given atomic center is slightly less than the nearest 
neighbor distance. In the present calculations, the use of overlapping spheres is 
essential, particularly for the small interatomic distances (r < 2AU) in the pseudo- 
molecules formed at the projectile energies used in the experiments. The projectile 
energy range of interest in Ne-Mg (1-3 KeV) corresponds to interatomic distances in 
the range 0.47 AU to about 1.0 AU (ref. 6). 
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INTERACTIVE SCF-Xa-SW CALCULATIONS 


The computer programs necessary to carry out SCF-xa-SW calculations in an 
interactive manner have been in^lemented on the LeRC IBM 370 con^uter using the 
TSS (derating system. In order to facilitate use of the programs interactively from a 
terminal, and to minimize the amount of data to be entered from the terminal, a series 
of special PROCDEFS procedure definitions) were created (appendix I). The use of 
these are demonstrated below, A con^lete xa calculation entails the following 
sequence of steps, each of which involves the use of a different program package: 

(a) Generation of atomic potentials for each atom or ion in the molecule using a 

modified Herman Skillman SCF program 

(b) Generation of a molecular potential by superposition of atomic potentials 

(c) Generation of a contracted molecular potential for use in steps (e) and (f) 

(d) Generation of symmetry input for use in steps (e) and (f) 

(e) Calculation of the energy levels of the molecule 

(f) SCF calculation starting with the results of step (e) plus subsequent interation 
In addition, one may calculate the wave functions and charge densities corresponding to 
the SCF energy levels and also plot these quantities. (Thus far, these programs have 
not been implemented. ) The output from each successive step is stored in a specified 
element on disk. The programs for each successive step are run under the control 

of an appropriate PROCDEF. The input data for a given step is automatically con- 
structed by the PROCDEF by combining the previously stored elements with minimal 
input entered from the terminal. 

To illustrate the use of the PROCDEFS referred to above, we shall consider the 
sequence of operations involved in carrying out an SCF calculation for a Mg-Ne 
pseudomolecule. 

A. Generation of atomic potentials 

The format of the input data cards necessary to generate atomic potentials using 
the modified Herman Skillman program are described in the unpublished operating 
instructions for the xa programs. A listing of the input cards for Ne, stored in 
element NEI, is given in appendix 11(B). To generate the atomic potential for Ne and 
to store the card image output in element NEI, one enters 

_HSRUNNEI, NEI 

at the terminal. The terminal will respond with 

DDNAME FT06F001 UNDEFINED ENTER Y. . . . 
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To print the condensed tabular output at the terminal, one enters D. The eigenvalues 
for each iteration, and the total charge density and core charge density for the final 
iteration are then printed at the terminal. Card image output is stored in element NEl. 

B. Generation of molecular potential 

The input data cards for a Ne-Mg pseudomolecule, stored in element NEMG, are 
listed in appendix 11(C). The cards are for the case where the Ne atom has coordinates 
(2, 0, 0) AU and Mg (-2, 0, 0) AU. To run the program, one enters 

_MRUNNEMG, VNM 

at the terminal. The PROCDEF loads NEMG and goes into the edit mode now assum- 
ing that the Ne and Mg total and core charge densities are stored in Nel and MGl, 
respectively, one next types in 

B; LOA NEl; B; LOA MGl 

This loads NEl and MGl. To reenter the PROCDEF and to continue the calculation, 
one types 


DEFAULT SYSINX=E 

The output molecular potential is stored in element VNM. For this program tabular 
output is listed on the high-speed printer, rather than at the terminal, because the out- 
put is quite lengthy (2200 lines for NeMg). 

C. Greneration of contracted molecular potential 

The next step is the generation of a contracted molecular potential assuming speci- 
fic values of radii for the atomic and outer spheres in the molecule. If the input mole- 
cular potential is stored in element VNM, the contracted potential is to be stored in 
VVM and the radii of the outer, Ne, Mg spheres are 0.705, 0.4332, 0.470 AU, 
respectively, one enters at the terminal: 

_CRUNT VNM, VVM 

I .705 .4332 .47 0.0 

DEFAULT SYSINX = E 


The terminal responds with 

DDNAME FT06F001 UIU)EFINED. . . . 
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upon entry of D, the calculation proceeds and tabular output is printed at the terminal. 
The table of the output contracted molecular potential is listed in appendix 11(E). 

D. Generation of symmetry input 

This program calculates spherical harmonic basis functions for a given set of 1 
values, and writes out card images which are stored in a specified element. The cal- 
culation must be repeated for each unique group of equivalent atoms and each symmetry 
group. The output elements must then be edited to form the input for E and F. Sup- 
pose the required input cards, described in the xa operating instructions, have been 
placed in element PIN and the output element name is PIO. One would enter 

_YLRUN PIN, PIO 

at the terminal, which then responds with 

DDNAME FT06001 UNDEFINED. . . . 

Upon entry of D, the calculation proceeds and tabular output is printed at the terminal. 
The final edited versions of the symmetry ouiput for the ct, tt representations of the 
group of the heteronuclear diatomic molecule are given in appendix II (A). 

E. Energy level calculation 

Suppose that one wishes to search for energy levels in the energy range -0. 5 AU 
to -10. AU in steps of 0. 5AU. One would enter 

_FRUN VNM, PIZ 

at the terminal. (The contracted potential is stored in element VNM, and the symmetry 
ii^ut is stored in PIZ.) The terminal will respond with 

DDNAME FT06F001 UNDEFINED. . . . 

and one enters D to print tabular output at the terminal. When the terminal indicated 
the editor is entered, one types in 

I 3, 1,. 5, -.5,-10., 0 

DEFAULT SYSINX=E 

and the calculation proceeds. The determinant value for each trial energy in the range 
specified is then calculated and printed. If the sign of the determinant changes in a 
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given energy range, the program automatically interpolates for the energy correspond- 
ing to zero determinant (that is, energy eigenvalue). 

F. SCF calculation 

The final step in the xoi calculational procedure is the iterative SCF calculation. 
This is performed in two stages under control of two different PRCKJDEFs. In the first 
stage, the necessary input data for the second stage is constructed automatically by the 
PROCDEF. This is convenient because of the complexity of the input data and because 
similar input data is needed for the SCF calculations at different interatomic distancea 
The PROCDEF DIS constructs SCF iiput for NeMg in the case where there are 4 
different sigma levels and 2 different pi levels occupied in the pseudomolecule. One 
must also specify the number of iterations to be performed and the energies of the 
frozen core levels. Suppose the number of iterations desired is 8, frozen core levels 
are found at -60.9204, -92.6720 AU, sigma states at -6.30841, -4.51048, -2.88924 
-1.46411 AU, and pi states at -3.63858 and -1.601967. One enters at the terminal 

DIGS, -60.9204, -92.6720, -6.30841, -4.51048, 

-2.88924, -1.46411, -3.63858, -1.601967 

This generates the SCF input, placing it in element DIG. The SCF input generated 
for this case is listed in appendix 11(F). Next, with input potential in VM and output 
in VMl, the SCF step is started by typing 

_SCFR VM, DIG, VMl 


The terminal responds with 

DDNAME FT06F001 UNDEFINED 

and one enters D to proceed. In each iteration the input levels and occupation numbers^ 
the various components of potential energy, the total potential energy, kinetic energy, 
total energy and virial coefficient are printed. Also printed are the total charge in each 
region, the largest difference between the old and new potentials, calculated energy 
eigenvalues for the old potential, and predicted eigenvalues for the new potential. After 
the last iteration, charges are given as fractions of the total charge within the atomic 
spheres, and decomposed in their spherical-harmonic (s, p, d, f, etc.) con^onents. 
Finally, the new potential is printed at the terminal. 
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APPENDIX I - PROCDEFS FOR RUNNING THE Xa PROGRAMS 
A. HERMAN SKILLMAN ATOMIC POTENTIAL PROCDEF 


LIST OtIOOO 

HSRUN 0000000 PROCDEF HSRUN 

HSRUN 0000100 PAR AM TN.OUT 

HSRUN 0000150 DEFAULT SYSINX=E 

HSRUN 0000200 IN 

HSRUN 0000300 7.E 

HSRUN 0000400 RED IT TN 

HSRUN 0000500 B?LOA INPLUVS INPU;0 

HSRUN 0000600 HRUN HSKK,OUT 

HSRUN 0000700 ABEND 


PROCDEF HRUN 
LIST 100.2500 

HRUN 0000100 PAR AM JBLIB.CARD 

HRUN 0000200 PPL I 

HRUN 0000300 DDEF ASD, VP, JBLIB, OPTION=JOBLIB 

HRUN 0000400 F I LEDEF FT05F00 1 , VS , I NPU , DCB= ( RECFM=F , LRECL=SO ) 

HRUN 0000500 FI LEDEF FT07F001 , VS, CARD, DCB= ( RECFM=F, LRECL=80 ) 

HRUN 0000600 OSRUN MAIN 

HRUN 0000700 PRINT CARD 

HRUN 0000800 FILEREL FT05F001 , FT07F001 

HRUN 0000900 ERASE I NPU 

HRUN 0001000 RELEASE ASD 


B. MOLECULAR POTENTIAL PROCDEF 


PROCDEF MRUN 
LIST 100,2500 

MRUN 0000100 PARAM TNP.OUT 

MRUN 0000200 DEFAULT SYSINX=E 

MRUN 0000300 IN 

MRUN 0000400 7.E 

MRUN 0000500 RED IT TNP 

MRUN 0000600 PAUSE 

MRUN 0000700 DEFAULT SYSINX=G 

MRUN 0000800 CONED IT 

MRUN 0000900 B;LOA I NPU; VS I NPU 

MRUN 0001000 Q 

MRUN 0001100 SRUN MOLE, OUT 

MRUN 0001200 ABEND 
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PROCDEF SRLIN 
LIST 100,2500 

SPUN 0000100 PARAM JLIB,CARD 

SPUN 0000200 PPL I 

SPUN 0000300 F I Lt'DEF FT06F00 1 , VS , OUT , DCB= < RECFM=F , LRECL® 136), RET=T 

SRLIN 0000400 HRUN JLIB,CARD 

SRUN 0000500 PRINT OUT 

SPUN 0000600 FILEREL FT06F001 


PROCDEF HRUN 
LIST 100,2500 

HRUN 0000100 PARAM JBLIB,CARD 

HRUN 0000200 PPL I 

HRUN 0000300 DDEF ASD, VP, JBLIB, OPTION=JOBLIB 

HRUN 0000400 FILEDEF FT05F001 , VS, INPU, DCB= ( RECFM=F, LRECL=80 ) 

HRUN 0000500 FILEDEF FT07F001 , VS, CARD, DCB= ( RECFM=F, LRECL=80 ) 

HRUN 0000600 OSRUN MAIN 

HRUN 0000700 PRINT CARD 

HRUN 0000800 FILEREL FT05F001 , FT07F001 

HRUN 0000900 ERASE INPU 

HRUN 0001000 RELEASE ASD 


C. CONTRAC PROCDEF 


PROCDEF CRUNT 
LIST 100,2500 

CRUNT 0000100 PARAM TNP^,OUT 

CRUNT 0000200 DEFAULT SYSINX=E 

CRUNT 0000300 IN 

CRUNT 0000400 7.E 

CRUNT 0000500 RED IT TNP-' 

CRUNT 0000600 T;LjA RNEMG;T;N 2;D 

CRUNT 0000700 PAUSE 

CRUNT 0000800 DEFAULT SYSINX=G 

CRUNT 0000900 CONED IT 

CRUNT 0001000 B;LOA INPU? VS INPU;Q 

CRUNT 0001100 HRUN CONT , OUT 

CRUNT 0001200 RED IT OUT 

CRUNT 0001300 T;P 1000?Q 

CRUNT 0001400 ABEND 
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PROCDEF HRUN 
LIST 100,2500 


HRUN 0000100 PAR AM JBLIB,CARD 

HRUN 0000200 PPL I 

HRUN 0000300 DDEF ASD, VP, UBLIB, OPTION=JOBLIB 

HRUN 0000400 FILEDEF FT05F001 , VS, INPU, DCB= (RECFM=F, LRECL=80) 

HRUN 0000500 FILEDEF FT07F001 , VS, CARD, DCB= (RECFM=F» LRECL«80) 

HRUN 0000600 OSRUN MAIN 


HRUN 0000700 PRINT CARD 

HRUN 0000800 FILEREL FT05F001 , FT07F001 

HRUN 0000900 ERASE INPU 

HRUN 0001000 RELEASE ASD 


D. SYMMETRIZATION PRCX5RAM PROCDEF 


PROCDEF YLRUN 


LIST 100, 1500 

YLRUN 0000100 PAR AM YLIN,OUT 

YLRUN 0000200 DEFAULT SYSINX=E 


YLRUN 0000210 IN 

YLRUN 0000220 7.E 

YLRUN 0000300 RED IT YLIN 


YLRUN 0000400 B;LOA INPU; VS INPU;Q 

YLRUN 0000500 PPL I 

YLRUN 0000600 TTPE FT08F001,SCR 


YLRUN 0000700 HRUN HAR,OUT 

YLRUN 0000800 FILEREL FT08F001 

YLRUN 0000900 RED IT OUT 


YLRUN 0001000 T;P 300; Q 

YLRUN 0001100 ABEND 


PROCDEF HRUN 


LIST 

100, 1900 

HRUN 

0000 1 00 

HRUN 

0000200 

HRUN 

0000300 

HRUN 

0000400 

HRUN 

0000500 

HRUN 

0000600 

HRUN 

0000700 

HRUN 

0000800 

HRUN 

0000900 

HRUN 

0001000 


PAR AM JBLIB,CARD 
PPL I 

DDEF ASD , VP , JBL I B , OPT I 0N= JOBL I B 

F I LEDEF FT05F00 1 , VS , I NPU , DCB= ( RECFM=F , LRECL=80 ) 
F I LEDEF FT07F00 1 , VS , CARD , DCB= ( RECFM=F , LRECL=80 ) 
OSRUN MAIN 
PRINT CARD 

F I LEREL FT05F00 1 , FT07F00 1 
ERASE INPU 
RELEASE ASD 


PROCDEF TTPE 
LIST 100,1300 

TTPE 0000100 PARAM TAPE.SCRAT 

TTPE OOO02OO FILEDEF TAPE, VS, SCRAT, DCB= ( RECFM=V, LRECL=S00 ) , OSRECFM=V 


S,RET=T 
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E. ENERGY CALCULATION PROCDEF 


PROCDEF FRUN 
LIST 100,1900 

FRUN 0000100 PARAM NAME,DIA 

FRUN 0000200 DEFAULT SYSINX=E 

FRUN 0000300 IN 

FRUN 0000400 V.E 

FRUN 0000500 RED IT NAME 

FRUN 0000550 T; I ENERGY 

FRUN 0000600 B?LOA DIA;B 

FRUN 0000700 PAUSE 

FRUN 0000800 DEFAULT SYSINX=G 

FRUN 0000900 CONED IT 

FRUN 0001000 I 0,0, 0,0, 0,0 

FRUN 0001100 B?LOA INPLUVS INPU;Q 

FRUN 0001200 YRUN ENER,ENEI 

FRUN 0001300 ABEND 


PROCDEF YRUN 
LIST 100,1900 

YRUN 0000100 PARAM JBLIB,CARD 

YRUN 0000200 PPL I 

YRUN 0000300 DDEF ASD , VP , JBL I B , OPT I 0N= JOBL I B 

YRUN 0000400 DDEF BSD, VP, NJLIB, OPTION=JOBLIB 

YRUN 0000450 DDEF CSD, VP, TMA, OPTION=JOBLIB 

YRUN 0000500 FILEDEF FT05F001 , VS, INPU, DCB= ( RECFM=F, LRECL=80 ) 

YRUN 0000600 F I LEDEF FT07F00 1 , VS , CARD , DCB= ( RECFM=F , LRECL=80 ) 

YRUN 0000700 OSRUN MAIN 

YRUN 0000800 PRINT CARD 

YRUN 0000900 FILEREL FT05F001 , FT07FOO 1 

YRUN 0001000 ERASE INPU 

YRUN 0001100 RELEASE ASD 

YRUN 0001200 RELEASE BSD 

YRUN 0001300 RELEASE CSD 


F. PROCDEF FOR GENERATING SCF PROGRAM INPUT 


PROCDEF DIS 
LIST 0,1000 

DIS 0000000 PROCDEF DIS 

DIS 0000100 PARAM NUM, N1 , Ml , A1 , B1 , Cl , D1 , El , FI 

DIS 0000200 DEFAULT SYSINX=E 

DIS 0000250 ERASE RIG 

DIS 0000300 REDIT DIG 

DIS 0000400 PO 200IC /8/NUM/;N;DU 12/Nl / ; P? N? DU 12/Ml /? P; N; DU 11/Al/; 

p;PO 2500; DU 1 i/bi/;p;n;du ii/ci/;p;n 
DIS 0000500 DU 1 1 /Di/ ; p; N; DU ii/ei/;p;b;du ii/fi/;p;fil rig;q 
DIS 0000600 abend 
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G. SELF-CONSISTENT FIELD CALCULATION 


PROCDEF 

SCFR 


LIST 100,3000 


SCFR 

0000100 

PARAM NAME,DIA,SCFOT 

SCFR 

0000200 

DEFAULT SYSINX=E 

SCFR 

0000250 

ERASE I NPU 

SCFR 

0000300 

IN 

SCFR 

0000350 

SCF 

SCFR 

0000400 

y.E 

SCFR 

0000500 

REDIT DIA 

SCFR 

0000600 

po ioo;loa name;n;d?b 

SCFR 

0001000 

I 0,0, 0,0, 0,0 

SCFR 

0001100 

t;loa I NPU; vs inpu;q 

SCFR 

0001200 

XRUN SCFOT 

SCFR 

0001400 

REDIT SCFOT 

SCFR 

0001450 

t; 1 1 1 

SCFR 

0001500 

T;p 1000; FIL SCFOT 

SCFR 

0001600 

PAUSE 

SCFR 

0001700 

DEFAULT SYSINX=G 

SCFR 

0001750 

CONED IT 

SCFR 

0001800 

FIL SCFOT ;Q 

SCFR 

0001900 

ABEND 


PROCDEF 

XRUN 

LIST 100,2500 

XRUN 

0000100 

XRUN 

0000200 

XRUN 

0000300 

XRUN 

0000310 

XRUN 

0000320 

XRUN 

0000355 

XRUN 

0000360 

XRUN 

0000365 

XRUN 

0000370 

XRUN 

0000375 

XRUN 

0000380 

XRUN 

0000500 

XRUN 

0000600 

XRUN 

0000700 

XRUN 

0000800 

XRUN 

0001000 

XRUN 

0001100 

XRUN 

0001200 

XRUN 

0001300 

XRUN 

0001400 


PARAM CARD 
PPL I 

DDEF ASD , VP , ENER , OPT I ON=JOBL I B 
DDEF BSD , VP , N JL I B , OPT I ON= JOBL I B 
DDEF DSD , VP , TMA , OPT I ON=JOBL I B 
DDEF CSD, VP,SCFL,OPTION=JOBLIB 
TTPE FT08F001,SC8 
TTPE FTlOFOOl.SCO 
TTPE FT11F001,SC1 
TTPE FT12F001,SC2 
TTPE FT13F001,SC3 

F I LEDEF FT05F00 1 , VS , I NPU , DCB= ( RECFM=F , LRECL=80 ) 

FILEDEF FT07F001 , VS, CARD, DCB= ( RECFM=F, LRECL=80) 

OSRUN MAIN 

FILEREL FT 

ERASE I NPU 

RELEASE ASD 

RELEASE BSD 

RELEASE CSD 

RELEASE DSD 


EROCDEF TTPE 
LIST 100,1200 

TTPE 0000100 PARAM TAPE, SCR AT 

TTPE 0000200 FILEDEF TAPE, VS, SCRAT, DCB= ( RECFM=V, LRECL=800 ) , OSRECFM=VS, RET==T 
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APPENDIX n - INPUT DATA FOR THE Xa PROGRAMS 
A. DIATOMIC MOLECULE SYMMETRY INPUT FOR 
ENERGY, SCF CALCULATIONS 


non piz 

L0MIIN6 PIZ 
EDIT 


r;p ioo;o 

TOP RECORD 
0000020 10 2 
0000200 1 1 
0000300 0 1 
000040C 2 1 
0000500 1 1 
0000600 3 1 
0000700 0 1 
0000800 3 1 

loOlW 4 1 
0001100 1 1 
0001200 4 1 
0001300 3 1 
0001500 1 1 
0001600 0 1 
0001700 2 1 
OOOIBOG 1 1 
0002000 1 1 
0002100 0 1 
0002200 2 1 
0002300 1 1 
EOF 


PI Z REP 
1 1.00000000 
1 1.00000000 
1 1.00000000 
1 1.00000000 
1 1.00000000 
1 1.00000000 
2 1.00000000 
2 1.00000000 
3 1.00000000 
3 1.00000000 


REPIT SI6HB 
LOADING S16HA 
EDIT 

Tip lOOiQ 

TOP RECORD 
0000020 B 2 

0000200 0 1 

0000300 0 1 

0000400 1 1 

0000500 1 1 

0000600 2 1 

0000700 0 1 

0000800 2 1 

0000900 2 1 

0002100 0 1 

0002200 0 1 

0002300 1 1 

0002400 1 1 

0003000 0 1 

0003100 0 1 

0003200 1 1 

0003300 1 . 1 

EOF 


SIGHA REP 

1 i.ooooooo: 
1 1.00000000 
1 1.00000000 
1 1.00000000 

2 1.00000000 

2 1.00000000 
3 1.00000000 
3 1.00000000 
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B. NEI; HERMAN SKILLMAN PROGRAM NE INPUT 


0000100 

0000120 

0000200 

CONT PUHC 

NEUTRAL NEON ATOMIC POTENTIAL 
0 .001 .00001 991 1 50 1 

; IS FROZEN 

.99128 

.93266 

.73081 

0000300 

1.00 

.99205 

.98387 

.97552 .96705 

.95850 

.99990 

.92907 

0000^00 

.91553 

.89862 

.88202 

.86578 .89999 

.83951 

.81997 

.80989 

.79058 

.77670 

0000500 

.76316 

.73706 

.71215 

.68835 .66558 

.69381 

.62303 

.60322 

.58938 

. 56698 

0000600 

.54948 

.51795 

.48923 

.96280 .93822 

.91523 

.39363 

.37332 

.35922 

.33627 

0000700 

.31941 

.28880 

.26198 

.23853 .21805 

.20015 

.18997 

.17068 

. 15851 

.19773 

oooosoo 

.13812 

.12180 

.10847 

.1 .1 

. 1 

. 1 

.1 

• 1 

• 1 

0000900 

.1 

.1 

.1 

.1 .1 

.1 

.1 

. 1 

• 1 

• 1 
1 

0001000 

.1 

.1 

.1 

.1 .1 

• 1 

• 1 

• 1 

• 1 

« 1 
1 

0001100 

.1 

.1 

.1 

.1 .1 

. 1 

. 1 

.1 

• 1 

• 1 

0001200 

.1 

.1 

.1 

.1 .1 

.1 

• 1 

* 1 

• 1 

• 1 
1 

0001300 

0001400 

0001500 

0001600 

0001700 

OOOISOO 

.1 

10 1 
100 2. 

200 2. 

210 6. 

.1 

2 

.1 

-60.92 

-3.168 

-1.471 

.1 .1 

.1 

.1 

• 1 

• 1 

• 1 




c. 

NEMG: INPUT 

0000100 

1 




0000200 

3 

3 

1 PCH 

PRT 

0000300 

OUT 

0 

2 1 

0. 

0000400 

2 

1 

2 3 

1 

0000500 

NE 

0 

1 2 

2. 

0000600 

3 

1 

3 


0000700 

MG 

0 

1 3 

-2. 

0000800 

2 

1 

2 



FOR MOLEC PROGRAM (NE-MG) 


1 .73 



0. 

0. 

.73 

3 

0. 

0. 

.73081 

0. 

0. 

.72913 


D. RNEMG: CONTRAC INPUT FOR NE-MG 


0000100 2 
0000120 4.946 

2 

2.19992 

2.946 

0000300 0. 





0000^00 0 . 
0000500 0. 


0 * 


E. VVM: CONTRACTED MOLECULAR POTENTIAL 
aNPUT FOR ENERGY, SCF PROGRAMS) 


1 1 - 0.0 0.0 

3312 0.73000 1 2 

OUT 0 1 0 93 4 0.0 0.0 

22 92 62 82 102 122 192 162 

0.5893213E+00 0 . 6197398E+00 0 .6551989E+00 
0.7613893E-t-00 0 .7968029E-^00 0 .8322166E->^00 
0.9389579E+00 0 . 9738711E+00 0 . 1009289E+01 
0.1115520E+01 0.1150931E+01 0 . 1186393E+01 
0.1292579E+01 0 . 1327991E+01 0. 1398816E+01 
0.1611293E+01 0.1682118E+01 0 . 1752999E+01 
0.1965921E+01 0 .2036296E+01 0.2107072E+01 


0.0 0.70500 0.73000 
82 202 

0.6905621E+00 0 .7259757E+00 
0.8676302E+00 0 . 9030938E+00 
0.1099696E-»-01 0 . lOSOlOSE-i-Ol 
0.1221755E+01 0 . 1257167E+01 
0.1969692E+01 0 . 1590967E+01 
0.1823770E+01 0 . 1899595E+01 
0.2177897E+01 0 .2298723E+01 
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0.23195<i9E+01 0 . 2390374E+01 0 . 296120 OE+01 0 . 2532025E+0 1 0 . 2602851E+01 11 

0.2673676E+01 0 . 2799502E+0 1 0 . 2886 155E+0 1 0 . 3027808E+01 0 . 3169961E+01 12 

0.3311119E+01 0 .3952767E+01 0 . 3599920E+01 0 . 3736 07 3E+0 1 0 . 3877727E+01 13 

0 .9019380E+01 0 .9161033E+01 0 .9302686E+01 0 . 9999339E+01 0 . 9585992E+01 19 

0 .9727695E+01 0 . 9869298E+01 0 .5010951E+01 0 .5152609E+01 0 . 5299257E+01 15 

0 .5935910E+01 0 . 5577563E+01 0 . 5860871E+01 0 . 619917 9E+01 0 . 6927987E+01 16 

0.6710795E+01 0 . 6999103E+01 0 . 7277911E+01 0 . 75607 1 9E+01 0 , 7899028E+01 17 

0.8127336E+01 0 .8910699E+01 0 . 86 93952E+0 1 0 .8977260E+01 0 . 9260568E+01 18 

0. 9593876E+01 0 . 9827189E+01 0 . 1 0 1 1 09 9E+02 0 . 1 0 39380 E+ 0 2 0 . 1067711E+02 19 

0.1096092E+02 0 , 1129372E+02 0 . 1 181 039E+02 0 . 1237696E+02 0 . 1299358E+02 20 

0.1351020E+02 0 . 1907681E+02 0 . 1969393E+02 0 . 1521005E+02 0 . 1577667E+02 21 

0.1639326E+02 0 . 16 90 985E+02 0 . 1797699E+02 22 

-0. 2963831 E+02-0. 2177 068E+ 02-0. 1933852E+02-0. 17257 06 E+02-0.1596275E+02 23 

-0. 1 3906 90 E+02-0. 12551 36 E+02-0.1136551E+02-0.1032933E+02-0.9907099E+01 29 

-0.8596968E+01-0 . 787785 lE+0 1-0. 7238875E+0 1-0. 666 9029E+0 1-0. 6159399E+01 25 

-0.5702175E+01-0.5290991E+01-0.9920012E+01-0.9589592E+01-0.9280350E+01 26 

-0. 9003827 E+01-0. 37 51899E+0 1-0. 3311 037 E+0 1-0. 2990 17 lE+0 1-0. 2625712E+01 27 

-0. 235721 3E+01-0. 21 26528E+0 1-0. 1927226E+01-0. 175917 OE+0 1-0. 16 03210E+01 28 

-0. 197 0959E+0 1-0. 13596 03E+0 1-0. 1251829E+01-0. 116 0689E+0 1-0. 107 9531 E+01 29 

-0.1006989E+01-0.9918797E+00-0.8832099E+00-0.8301987E+00-0.7819797E+00 30 

-0.7380802E+00-0.6979955E+00-0.6272957E+00-0.5670267E+00-0.5151798E+00 31 

-0. 97 0 0 997E+00-0. 930567 1 E+0 0-0. 3956957E+00-0. 3695869E+0 0-0. 3367 962E+00 32 

-0.3117971E+00-0.2892009E+00-0.2686865E+00-0.2999905E+00-0.2328922E+00 33 

-0. 2172096 E+0 0-0 .20277 18E+0 0-0. 1899593E+00-0. 1771 599E+0 0-0. 1657573E+00 39 

-0.1551822E+00-0. 1953576E+00-0.1276935E+00-0.1123301E+00-0.9892082E-01 35 

-0. 87 18550E-0 1-0. 7689518 E-0 1-0. 6785882E-01-0. 599121 lE-0 1-0. 5291689E-01 36 

-0.9675392E-01-0.9132001E-01-0.3653039E-01-0.3230766E-01-0.2858673E-01 37 

-0. 2529121 E-0 1-0. 2237 992E-0 1-0. 1 980313E-01-0. 1755022E-0 1-0. 1559865E-01 38 

-0. 1 377395E-0 1-0. 1221991 E-0 1-0. 9601 095E-02-0. 7598936 E-02-0.5935065E-02 39 

-0.96 7567 9E-02-0. 37 05239 E-0 2-0. 2968737 E-02-0. 238 9 991 E-02-0. 186 2353E-02 90 

-0.1335975E-02-0.6751183E-03-0.3205158E-03 91 

NE 010 1 72 69 0.90000 0.0 0.0 0.93320 0.73081 92 

20 90 60 160 160 160 160 160 160 160 93 

0.2059696E-02 0 . 91 09390E-02 0 . 6 1 69 081 E-02 0 . 8218773E-02 0 . 1027396E-01 99 

0.1232816E-01 0 . 1938285E-01 0 . 1693759E-01 0 . 189 9223E-0 1 0 . 2059692E-01 95 

0.2260161E-01 0 .2965630E-01 0 . 2671099E-01 0 .2876569E-01 0 . 3082038E-01 96 

0.3287507E-01 0 . 3992976 E-0 1 0 . 36 98995E-0 1 0 .3903919E-01 0 . 9109383E-01 97 

0.9520322E-01 0. 993126 1 E-0 1 0 . 5392200E-01 0 .5753139E-01 0 . 6169078E-01 98 

0.6575012E-01 0 . 6985950E-01 0 . 7396889E-01 0 .7807827E-01 0 .8218765E-01 99 

0.8629709E-01 0 . 9090692E-01 0 . 9951580E-01 0 . 9862518E-01 0 . 1027396E+00 50 

0.1068939E+00 0 . 1 1 0 9533E+0 0 0 . 1150627E+00 0 . 1191721E+00 0 . 1232815E+00 51 

0.1315002E+00 0 . 1397190E+00 0 . 197 9378E+00 0 . 1561565E+00 0 . 1693753E+00 52 

0.1725991E+00 0 . 1808128E+00 0 . 1890316E+00 0 . 1972509E+00 0 . 2059691E+00 53 

0 .2136879E+00 0 . 2219067E+00 0 . 2301259E+00 0 . 2383992E+00 0 . 2965630E+00 59 

0 .2597817E+00 0 . 2630 0 05E+0 0 0 . 2712193E+00 0 .2799380E+00 0 . 2876568E+00 55 

0.3090993E+00 0 . 3205318E+00 0 . 3369699E+00 0 .3539069E+00 0 . 3698999E+00 56 

0.3862820E+00 0 . 9027 1 95E+ 0 0 0 . 9191570E+00 0 .9355996E+00 0 . 9520321E+00 57 

0.9689696E+00 0 . 989 9 07 2E+ 0 0 58 

-0 . 96996 91 E+09-0 .982757 9E+09-0 . 3205081 E+09-0 . 23937 60 E+0 9-0 . 1906919E+09 59 

-0 .1582328 E+ 09-0 . 1350958 E+09-0 . 1176599E+09-0 . 109 1272E+09-0 . 9330597E+03 60 

-0. 8995159 E+0 3- 0. 77 0 7 38 3 E+0 3-0. 7 083189E+ 03- 0.659829 5 E+0 3- 0.6 089727 E+0 3 61 

-0.5679255E+03-0 .532159 9E+0 3-0. 50 037 99E+03-0. 97 19570 E+0 3-0. 996 3889 E+0 3 62 

-0 .9022603E+03-0 . 36 55332 E+0 3-0 . 3 395020 E+03-0 . 307 9980 E+0 3-0 . 289 977 3E+03 63 

-0. 269 9182E+0 3-0 .297 257 9 E+0 3-0. 231 596 3E+ 03-0. 21 7 6 178 E+0 3-0 .20 506 96 E+03 69 

-0. 1937971 E+03-0. 1839826 E+0 3-0. 1791375E+03-0.1655969E+03-0.1577623E+03 65 

-0.1505529E+0 3-0. 1938982E+03-0. 1377 389E+03-0. 1320216 E+0 3-0. 1267 030E+03 66 

-0. 11710 92E+03-0.1086993E+ 03-0 . 1012729E+03-0 . 996 7 226 E+02-0 . 8877 127 E+02 67 

-0 .8396828E+02-0 . 7868033E+02-0 . 793391 OE+02-0 . 7038782E+02-0 . 6677 909E+02 68 

-0 .6397263E+02-0 .6093959E+02-0.5763569E+02-0.5505083E+02-0.5265831E+02 69 

-0.5093915E+02-0. 9837 672E+02-0 .9695638E+02-0 .996 651 6 E+02-0 . 9299158E+02 70 

-0.39 95799 E+02-0. 37 289 36 E+0 2-0. 39 916 37 E+02-0. 328081 9E+02-0. 30 9227 9E+02 71 

-0.2923012E+02-0 .2770522E+02-0.2632792E+02-0.2507995E+02-0.2399675E+02 72 

-0.2291695E+02-0.2197952E+02 73 

0.2975518E-01 0 . 1192313E+00 0 . 2966817 E+00 0 .9209120E+00 0 . 6312920E+00 79 

0.8729719E+00 0 . 1 1 39896 E+0 1 0 . 1929028E+0 1 0 . 1736062E+01 0 .2057357E+01 75 

0.2389699E+01 0 . 2729987 E+0 1 0 . 3075709E+01 0 . 3929397E+01 0 . 3773910E+01 76 

0.9122313E+01 0 . 9967883E+01 0 .9809093E+01 0 . 5199590 E+0 1 0 . 5973190E+01 77 

0.6105685E+01 0 . 6699875E+01 0 . 7250917E+01 0 . 7755533E+0 1 0 . 821 1730E+01 78 

0.8618589E+01 0 .8976037E+01 0 . 9289736E+01 0 . 9595879E+01 0 . 9761108E+01 79 

0.9932391E+01 0 . 1006199E+02 0 . 1015215E+02 0 . 1020552E+02 0 . 1022960E+02 80 

0.1021198E+02 0.1017021E+02 0 . 1010181E+02 0 . 10 0 0 925E+02 0 . 9899886E+01 81 

0.9609775E+01 0 . 9263920E+01 0 .8870938E+01 0 .8995575E+01 0 . 7998799E+01 82 

0.7590126E+01 0 .7077803E+01 0 . 6618909E+01 0 . 61 6727 0 E+0 1 0 . 5728602E+01 83 

0.5305622E+01 0 . 9 90 07 0 1 E+0 1 0 . 9515993E+01 0 . 9151 058E+0 1 0 . 3807956E+01 89 

0.3986393E+01 0 . 3186051E+01 0 . 2906652E+01 0 . 2697 525E+0 1 0 . 2907896E+01 85 

0.1983599E+01 0 . 1625896E+01 0 . 1326861E+01 0 . 107856 1 E+0 1 0 . 8736291E+00 86 

0.7053668E+00 0 . 5678670E+00 0 . 9'559672E+00 0 . 3652385E+0 0 0 . 291 9209E+00 87 

0.2328500E+00 0 . 1853889E+00 88 
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MG 012 2 76 73 -0.<i0000 0.0 0.0 0 .^7000 0.72913 89 

20 40 60 160 160 160 160 160 160 160 90 

0.1933543E-02 0 . 3867 085E-02 0 .5800623E-02 0 . 7734161E-02 0 . 9667698E-02 91 

0.1160124E-01 0.1353477E-01 0 . 1546831E-01 0 . 1740185E-01 0 . 1 933539E-01 92 

0.2126892E-01 0 . 2320246E-01 0.2513600E-01 0 . 2706954E-01 0 .2900307E-01 93 

0.3093661E-01 0 . 3287015E-01 0 . 3480368E-01 0 . 3673722E-01 0 .3867076E-01 94 

0.4253784E-01 0 . 4640492E-01 0.5027201E-01 0 . 5413909E-01 0 .5800617E-01 95 

0.6187325E-01 0 . 6574029E-01 0 . 696 0732E-01 0 . 7347435E-01 0 .7734138E-01 96 

0.8120841E-01 0.8507544E-01 0 .8894247E-01 0 . 9280950E-01 0 . 9667653E-01 97 

0.1005436E+00 0 . 1044106E+00 0 . 1082776E+00 0 . 1121446E+00 0 . 1160117E+00 98 

0.1237457E+00 0 . 1314798E+00 0 . 1392139E+00 0 . 1469479E+00 0 . 1546820E+00 99 

0.1624160E+00 0 . 1701501E+00 0 . 1778842E+00 0 . 1856182E+00 0 . 1933523E+00 100 

, 0.2010863E+00 0 . 2088204E+00 0 . 2165545E+0 0 0 . 2242885E+00 0 . 2320226E+00 101 

0 .2397566E+00 0 . 2474907E+00 0 . 2552248E+00 0 .2629588E+00 0 .2706929E+00 102 

0.2861611E+00 0.3016294E+00 0 . 3170976E+00 0 . 3325658E+00 0 .3480341E+00 103 

0.3635023E+00 0 . 3789706E+00 0 . 3944388E+00 0 .4099070E+00 0 . 4253753E+00 104 

0.4408435E+00 0 . 4563118E+00 0 . 4717800E+00 0 . 4872482E+0 0 0 .5027165E+00 105 

0.5181847E+00 106 

-0.12357 91 E+0 5- 0.6 151 395E+04-0. 408240 OE+04-0. 304781 OE+04-0. 2427 OOOE+04 107 

-0. 20 130 93 E+ 04-0. 17 17426E+04-0. 149566 9E+04-0. 1323 191 E+04-0.1185215E+04 108 

-0 . 1072334 E+0 4-0 . 97828 03E+03-0 .8987 10 9E+0 3-0 . 830 5254 E+03-0 . 77 14485E+03 109 

-0.7197749E+03-0.6741995E+03-0.6337068E+03-0.5974961E+03-0.5649253E+03 110 

-0. 5087217 E+03-0. 46 19558E+03-0.4224521E+03-0.3886562E+03-0.3594275E+03 111 

-0. 33390 92E+03-0. 31 14468E+03-0.2915303E+03-0.2737568E+03-0.2578042E+03 112 

-0. 24341 11 E+03-0. 2303642E+03-0. 2184866E+03-0 .2076313E+03-0.1976746E+03 113 

-0 . 18851 19 E+0 3-0 . 18 00540 E+03-0 . 1722247 E+03-0 . 16 4 9581 E+0 3-0 . 1581 972E+0 3 114 

-0. 1460 006 E+03-0. 1353075E+03-0.1258639E+03-0. 11746 95E+03-0.1099650E+03 115 

-0 . 1032213E+03-0 . 97 1 3354E+02-0 . 9161501 E+02-0 . 8659364E+02-0 . 8200900E+02 116 

-0.7780 9 97 E+0 2-0. 7 395287 E+02-0. 7040 030 E+02-0. 67 11993E+ 02-0. 6408366 E+02 117 

-0 .6 1266 98 E+0 2-0 .586 4844E+02-0 .5620 914E+02-0 .53 93243E+0 2-0 .518036 OE+02 118 

-0. 47 938 90 E+02-0. 44527 13E+ 02-0. 4149831 E+02-0. 387 96 22E+02-0. 3637531 E+02 119 

-0.3419826E+02-0.3223416E+02-0.3045724E+02-0.2884576E+02-0.2738127E+02 120 

-0 .26 047 94 E+02-0 .24832 18E+02-0. 2372220 E+02-0 .227 0770 E+02-0. 2177 972E+02 121 

-0.2093042E+02 122 

0.4580327E-01 0 . 1749096E+00 0 . 3757172E+00 0 . 6376935E+00 0 . 9512933E+00 123 

0.1307882E+01 0 . 1699668E+01 0 . 21 19631E+01 0 .2561462E+01 0 . 3019513E+01 124 

0 .3488731E+01 0 . 3 96 46 1 9E+0 1 0 . 4443184E+01 0 . 49208 97 E+0 1 0 . 5394655E+01 125 

0.5861739E+01 0 . 6319783E+01 0 . 6766746E+01 0 .7200876E+01 0 . 7620689E+01 126 

0.8412610E+01 0.9135054E+01 0 . 9783422E+0 1 0 . 1 035542E+02 0 . 1085061E+02 127 

0.1126999E+02 0 . 1161568E+02 0 . 1189066E+02 0 . 1209854E+02 0 . 1224336 E+02 128 

0 . 1232944E+02 0 . 1236 126E+02 0 . 1234337E+02 0 . 1228029E+02 0 . 1217645E+02 129 

0. 1203616E+02 0 . 1 186354E+02 0 . 1166253E+02 0 . 1143685E+02 0 . 1119000E+02 130 

0 . 1064555E+02 0 . 1 0 05243E+02 0 . 9430264E+01 0 . 87 95227 E+0 1 0 . 8160419E+01 131 

0 .7536180E+01 0 . 6 930467E+0 1 0 . 6349185E+01 0 .5796524E+01 0 . 5275245E+01 132 

0.4786945E+01 0 . 4332289E+01 0 . 391 1202E+0 1 0 . 3523043E+0 1 0 . 3166745E+01 133 

0.2840932E+01 0 . 25440 16E+0 1 0 .2274278E+01 0 .2029923E+01 0 . 1809141E+01 134 

0 . 1431144E+01 0 . 1126569E+01 0 .8829421E+00 0 .6893017E+00 0 .5362442E+00 135 

0 .4158553E+00 0 . 3215727E+00 0 .2480209E+00 0 . 1908406E+00 0 . 1465265E+00 136 

0 . 1122804E+00 0 . 8588272E-0 1 0 . 6558204E-01 0 . 50 0 0 32 9E-0 1 0 . 3807138E-01 137 

0.2894885E-01 138 

-0.1598589E+02 139 


F. RIG: SCF INPUT (TRIAL EIGENVALUES AND SYMMETRY CARDS) 


0000100 0 1 

0000200 6, 8, 1,1, 5, 6, 0,0, 2, .05,1. 

0000300 2. -60.9204 

0000400 2. -92.6720 

0000500 SI 2. -6.30841 .005 110 

0000600 8 2 SIGMA EP 

0000700 0 1 

0000800 oil 1.00000000 

0000900 1 1 

0001000 1 1 1 1.00000000 

0001100 2 1 

0001200 oil 1.00000000 

0001300 2 1 

0001400 211 1.00000000 

0001500 0 1 

0001600 012 1.00000000 

0001700 1 1 
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0001800 

1 


1 

2 

1.00000000 




0001900 

0 


1 






0002000 

0 


1 

3 

1.00000000 




0002100 

1 


1 






0002200 

1 


1 

3 

1.00000000 




0002300 

52 

2 . 

- 

4.51048 

.005 

1 

1 

0 

0002^00 

S 3 

2 . 

- 

2.88924 

.005 

1 

1 

0 

0002^20 

54 

2 . 

- 

1.46411 

.005 

1 

1 

0 

0002500 

PIl 

4 . 

- 

3.63858 

.005 

5 

1 

0 

0002600 

10 


2 

PI Z 

REP 




0002700 

1 


1 






0002800 

0 


1 

1 

1.00000000 




0002900 

2 


1 






0003000 

1 


1 

1 

1.00000000 




0003100 

3 


1 






0003200 

0 


1 

1 

1.00000000 




0003300 

3 


1 






0003^00 

2 


1 

1 

1.00000000 




0003500 

4 


1 






0003600 

1 


1 

1 

1.00000000 




*0003700 

4 


1 






0003800 

3 


1 

1 

1.00000000 




•0003900 

1 


1 






0004000 

0 


1 

2 

1.00000000 




*0004100 

2 


1 






.0004200 

1 


1 

2 

1.00000000 




0004300 

1 


1 






0004400 

0 


1 

3 

1.00000000 




0004500 

2 


1 






0004600 

1 


1 

3 

1.00000000 




0004700 

PI 2 

4 . 

- 

• 1.601967 .005 

5 

1 

0 
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data to be entered manually at a remote terminal. Listings of the special PROCDEFS and of 
typical input data are given. 
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